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In this repoirt a method is presented for the u 
calculation of the profile drag of airfoil sections, 

The nethod requires only a knowledge of the theoretical 
velocity diatr ihut ion and can he applied readily 'once 
this distribution is ascertained. Comparison of calcu- 
lated and experimental drag characteristics for several 
airfoils shows a satisfactory agreement. Sample calcu- 
lations are included. 



INTSODUOTION 



H. B. Squire and A. D, Young developed a method, 
presented in reference 1, for calculating the profile 
drag over an arbitrary airfoil. Ihis method requires 
a knowledge of the velocity distribution over the air— . 
foil, and of the estimated position of transition from 
laminar to turbulent flow. The velocity distribution 
can be obtained by the method of reference 2. Using 
these data the growth of the boundary layer can be com- 
puted, the region of the laminar flow being anajlyzed by 
means of a step— by— step solution of Pohlhausen" s equation, 
and the turbulent region by a similar method of solution 
of an equation developed in reference 1. from the thick- 
ness of the boundary layer and the local velocity at the 
airfoil trailing edge. Squire and Young evaluate the 
profile drag by means of a semi— empirical formula. {This 
procedure is very labcrious because of the step— by^step 
computations involved, and its application is limited 
by lack of a method for estimating the position of 
transition. 

£• Jacobs and A. X. von Soenhoff (reference 3) 

give formulas for estimating the transition position 
for velocity distributions having maximum velocity far 
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■back, and f. ^J!nPaFBraE"a!Rr^5^^RiJ.iiKan (reference 4) 
indicate a means for o'btaining this information for 
velocity distributions with maximum velocity- far forward. 
In these references methods are presented for directly 
calculating the houndary— layer growth in the laminar 
region^ These methods, taken in conjvtnction with the 
concise graphical method for analyzing the region of 
turbulent flow» which is developed in the present re~ 
port, make it possible to calculate quickly the profile-^ 
drag coefficient of an arbitrary airfoil. This proce- 
dure is outlined in a readily usable form. Several 
illustrative examples are worked in detail aad the cal- 
culated results are compared with experimental measure- 
ments. 



THBOaT 



Ihe profile drag of an airfoil in free air is ao— 
count-ed for by the momentum loss in the wake far down- 
stream. It is possible to 'calculate this momentum loss 
by tracing the boundary^layer growth along the airfoil 
surface and its development, as the wake, after leaving 
the airfoil trailing edge. Ihe growth of the boundary 
layer is a function of the chordwise velocity distribu-' 
tio» over the airfoil. In studying this growth the 
essential factors are: 

1. Position of transition from laminar to turbu- 
lent flow 

2» Boundary^layer thickness at transition 

3. Qrowth of the turbulent boundary layer 

Once the thickfless of the boundary layer at the trail- 
ing edge has been computed the pressure recovery In the 
wake can be accounted for by the seml—emplr leal theory 
of Sq^uire and Young. Shese va'Tious steps will now be 
considered in detail. 

1. Ira^sition point .- The profile drag of airfoils 
is markedly affected by the position of transition 
from laminar to turbulent flow, For smooth airfoils in 
low tnirbulence streams at Reynolds numbers of the order 
of 5,000,000, -transition generally occurs close to the 
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laminar separation point. Surface roughness, stream 
turljulence, or an extr.emely higb. Reynolds num'ber can 
move transition forward of this point, ^or airfoils 
having maximum velocity far forward ^ however, there is 
no practical method available for predicting the extent 
of this movement of transition ahead of the separation 
point. In such eases therefore, transition position 
will be assumed to coincide with separation. It is con- 
venient to make this simplification hecause reference 4 
contains a direct method for calculating the separation 

point of airfoils having fflaxioua velocity far forward, 
This :;iiethod approximates the forward portion of the 
chorduise velocity distr it)ut ion by two straight lines. 
She region ahead of maxintua velocity Is represented by 



(1) 



and the region directly behind maximum velocity by 



X = (a + h) 5a ~ h^-s.^ (2) 

»0 C V 0 / 



where the subscript m designates the coordinate of 
the naximum velocity point. It can be shown from refer- 
ence 4 that the vel.oclty at separation Tg is a func- 
tion of the slopes of the approximating straight lines 
and the maxirnxun velocity Vm; namely, 



where is the function of b/a plotted in figure 1. 

Por low-drag airfoils which have maximum velocity 
far back, separation occurs so closely behind maximum 
velocity that transition may be assumed to lie at that 
position, except at very high Reynolds numbers where it 
is possible for transition to move' ahead of the maximum 
velocity point. As is pointed" out in reference 3, exist- 
ing experimental data taken in flight indicate that 
transitib^ will oecur^in a region of increasing velocity 
if the boundary— layer 'Seynolds number reaches a value 
of the order of 8000 to i^OO. It Is possible to apply 
this orltorion by using the relation given in reference 3 
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for cstlcuj-ating the Tjoundary-layer Reynolds number in 
a region of increasing velocity, 

where 

Hg a — i. Tjoundary— layer Reynolds nunber 

S = — airfoil Reynolds number 

fi distance normal to the surface from the surface 
7 0 the point in the boundary layer at which 
the velocity is 0.707Vj 

velocity outside boundary layer at the point for 
which boundary layer is being computed 

Vq free— stream velocity 

y velocity outside the boundary layer at any value 
of x/c 

X distance along airfoil surface from leading edge 

c airfoil choi^rd 

By the use of equation (3) or (4) it is possible 
to find the approximate location of the transition 
point for most airfoil velocity distributions. 

2, Boundary-^ aver thipkness at transition .-- The 

* 

growth of the turbulent bound-ary layer depends upon the 
initial boundary— layer thickness after transition, 
iFlight experiments indicate that transition occurs 
abruptly at a point. ' Since there is no appreciable 
drag over such a small region, the momentum thlolcness 
of the boundary layer, 6, remains unchanged in pass- 
ing from laminar to turbulent flow. Horeovert for a 
Blasius velocity profile the momentum thickness is re- 
lated to the boundary^layer thickness by 

e = 0.389 6 (5) 
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Therefore if the thickness of the laminar boundary layer 
at transition is coaputed'," t'he initial value of the tur— 
.Ijulont nonentum thickness is also known. 

Sq^tiiro and Toung studied the momentxun thickness in 
terms of the factor, t, . The relation, given in refer- 
ence 1, 'between these quantities is 

r . -\ 

t = 2.56 logg I 4.075 I (6) 



The laminar "boundary— layer thickness at transition 
can be calculated by the methods of references 3 and 4, 
From reference 4 it can be shown that when the velocity 
distribution caJi be approximated by the straight lines 
of equations (l) and (2), then 



V 1 . 79 V c > y S 



where Ag, the function of the ratio of the slopes 

of the straight lines, b/a , is plotted in figure 1. 
Substituting equations (5) and (7) into equation (6) 
gives ne the approximate initial value of the i factor 
when the position of maximum velocity is far forward 

= 2.56 log^ 1^0.66 As (8) 

In reference 3 the relation given for the laminar 
boundary-'layer thickness in a region of increasing 
velocity is 

s Jo VoJ 

Therefore, for this case of transition occurring in a 
region of increasing velocity, equation (6) takes the 
form 

L = 2.56 log- Tl.lsY'^'^ ^Ji'^ e"! (10) 
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Equation (8) for the "double— roof " velocity dis- 
tributions and equation (lO). for "low-drag" velocity 
distributions are adequate for finding the initial value 
of the i factor for most airfoils. 

3. (Turbulent' boundary-layer growth .— The analysis 

of the turbulent boundary layer is based on equation 
(ll) (p. 8 of reference l). 



d 



where ' ? 

I = • — , the 5 factor 

To intensity of skin friction 

Prime indicates derivative with respect to (x/c) 

Squire and Young use a step— by-step xntegration to solve 
this equation, Ihis process is very laborious and can 
be simplified to pk direct graphical evaluation without 
loss of accuracy, The simplification is as follows: 

Divide the velocity distribution in the turbulent 

region at those positions where its slope changes rapidly. 
Approxiuate each of these segments by a curve of the form 



(V/Vo) 



where IZ and L are constants determined by making 
curve 12 pass through the end points of the velocity- 
distribution segment under consideration. For a region 
of decreasing velocity, substituting equation ( 12) in 
equatio.n (llj leads to 



- — V o ^ 
KPR+6.13 ~ Y ^ 



dS 
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where 



J = 10.41 t"^ e-°- = "*^ 



It is also possible to separate the varia^bles of 
equation (ll) "by replacing the actual velocity distri— 
"bution m each segment hy its average velocity and 
averag,e gradient. Ihis alternative method leads to an 
integral of form similar to the left— hand side of equa- 
tion (13) and gives the same value for the t, factor. 
There is, therefore, no need for further considering this 
alternative procedure, 

She left side of equation (13) cannot "he integrated 
"by elenentary means. By numerical integration 



Up 011+ 6.13 

has "been evaluated for the usual range of t, and ER. 
These values are given in table I and plotted in figure 2. 

Therefore, when Sj, Z, fl, and logg ^® known, the 

value of ts may he found from figure 2 "by using equa- 
tion (13) m the form 



I 



Equation (14) applies only in regions over which 
the velocity decreases, When the velocity remains con- 
stant, equation ( 11 ) simplifies to 
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H ff)^^^^ = -JL. f^' e^-"-£ (15) 
\.Vo/ c ^ 10.41 J^^ 

= © ^iO. 245468 - 1.254t a + 3.304 

-e ^I0.2455i - 1.254 ^1 + 3.204 

L J 

vrhich can 1)6 solved by the use of figure 3. 

In a region, over which the velocity is increasingt 
equation (14) is replaced by 

values for the integral appearing in equation (16) have 
been evaluated and are given in table II and figure 4. 

It is rarely necessary to divide the velocity dis— 
tribtition in the region of turbulent flow into any more 
than three segments. In each of these segments either 
equation (14), (15), or (l6) will apply. Therefore the 
incj?ease of J through the turbulent region can be 
calculated directly by the use of figures 2, 3, and 4. 

Calculation of the , drag coefficient .- Squire and 

Toung have shown in reference 1 that when the final 
value of the momentum thickness, Sji, and the relative 

velocity at the trailing edge of an airfoil are known 
the drag coefficient can be calculated by the relation 



where the subscript 37 indicates values at the trail- 
ing edg,e. The previous analysis has been m terms of 
the ^ factor; so, to avoid having to convert to terns 
of moncntuia thickness, it is possible to use equation 
(6) rewritten as 
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T8 0.39l4t , , 

= 0.3454 0 = \ (18) 

V 

Su'bstituting into eq[uation (17) gives as the profile- 
drag coefficient 

In figure 5, \ is plotted as a function of i. 
Therefore, when the value of i ,^ has teen computed, 

can "be read from figure 5, and the profile drag can then 
be calculated "by equation (19). 

APPlICAOJIOir 



The nethod for calculating the profile drag, which 
is outlined in the previous section, can best be illus- 
trated by indicating the actual steps in the calculation 
of several examples, 

Bsanple 1. Low>-drag airfoil at high lift .- The 

velocity distribution for the HACA 35—215 airfoil was 
computed at a lift coefficient of 1.23 by the method of 
reference 3 and Is given by the solid line in figure 6. 
The trailing-edge velocity was obtained by extrapolating 
the calculated velocity gradient near the trailing edge. 
Distances along the airfoil surface are nearly the same 
as distances along the chord except at the leading edge. 
Since at the leading edge the method of reference 4 for 
laminar boundary layers applies, and this method is ♦not 
significantly affected by replacing distances along the 
surface by their chordwise component, it is possible to 
make the simplification o:^ using chordwise distances in- 
stead of actual lengths along the airfoil surface. 3?his 
simplification is usually permissible, . 

The computed velocity distribution has stagnation 
point at x/c = 0„0270 on the lower surface and a maxi- 
mum velocity of 2.39S times free stream at x/c = 0.0037 
on the upper surface. Behind this point the velocity 
reaches Y/Vp = 3,210 at x/c = 0.0108. Therefore, 
from equation (l) 
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3.392 

* " 0.0270 + 0.0037 

and equation (2) 

■b - 8 '39 2 - 2.210 
0.0108 - 0.0037 

frorl figUBo-1, since "b/a = 0.33, 

$s = 1«054, = 6.0 

e<iuatioix (3) gives as the velocity at transition (separa- 
tion) 



= 78,0 



= 35.6 



which corresponds to a vtilue of x/c = 0,0090, If the 
Reynolds number is taken to be 5.75 x lO®, the i 
factor at transition by .eq[uation (S) is 

ti = 2.56 log© r 0.66 (6) (O.OZOf) JsS -/ET^ '^^'^1' ^^'^ 

The positions in the turbulent region of the upper sur- 
face at which the velocity gradient changes abruptly are 
appro::inately: x/c » 0.0727, V/Vq = 1.725 and 

x/o = 0,50, T/Tq = 1.433. Using 7/Vo « 0.768 as the 
trailing— edge velocity, it is now possible to trace the 
growth of f. By eciuation (12) 



IC,H = (Q.0727 - 0.0090) 5^75 ^ ^^f ^ g^g^ ^ ^^e 

\1.725> \2,27j 

and then from equation (14) and figare 2 



I, = 0.042 + lege • = 0.316 

* si. 725/ 



so that from figure 2 
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i , » 24,05 

x/o = 0.073 



Similarly 

2,3 = ( 0»500^- 0.073) ^ .^^ ^ ^ ^^^^ ^ ^^s 

\J7432y 



" ^ 1 N /' 1 N 

\1.735/ 

Ig = 0.060 + loge (^^^ = 0.347 



x/c = 0,50 

Z R (1.00 -^0,50^ 5^75 X 10« = 4.77 X 10^ 
\0.768/ Vl.4387 

I3 = 0.593 + loge ( ^'^^^^ » 1.216 

L '=31.90 ' • , 



Proa figure 5 . ' 

= 6.50 X 10* 

The lower surface has a maximxm velocity of ', 
Y/Yq - 0.963 at x/c = O.'SS. Assuming that transition 
occurs at the point of maximum veloQity, the initial 
t factor 18 found from eoLUation ('9) 



5.3 / 1 y-^' 



J ,aj6 ' (0.158) = 0.206 X 10 

5.75X10^ ^0.963, • 

and eouation (10) 

t, = 2.56 loge F 1.18 /oTioTTTo^ (0.963)5.75 X 10®"'= 20.5 

The final t, factor is found hy equattons (12), (14), 

and figure 2, 
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/ 1 ^ _ /"...X--^ 

\0.768/ S^O.963 / 

I = 0.022 + loe/5iii5^ .= 0.248 . 
\0.768/ 

and 

tj » 26.2 

from figure 5 

\fji = 0*69 X 10* 

rinally, sul)stituting the valu« of Xj of the upper 
plus the lower surface in equation (19) gives 

= 3(6.50 0.69) X 10* ^ 



Example 2. Transition ahead of maximum velooity . - 

As another illustration of the methods for calculating 
the profile drag, the drag of the experimental velocity 
distriMtion for the NAOA 66,2-430 airfoil at zero angle 
of attack, given in reference 5, will te calculated for 
transition fixed at x/c - 0.10 on both surfaces at a 
Reynolds numher of 6.35 X 10®, corresponding to a 
Mach nvimher of 0.194, The upper and lower surfaces are 
treated as "being divided into two turhulent regions each. 
These divisions are shown on .the plot of the experimental 
velocity distribution (fig. 7). The analysis is based 
on eo.uation (9), (lO), (12), and (16), lor the upper 
surface: 

\yj 6.35 (5.3) 
Si = 2,56 loge 1^1.18 (0.199X10"') 1.2 (6.35) 10*1= 19,2 

S>^ Y^°1^Y ^'^^^ ^^^^5^10^ = 41.2 X' 10^ 
S.1.20 )~ \,1.32/ 



I, = 0.001 + log ^i-'^^) = 0.095 
^ K, ,,1.30/ 



fron figure 4 



t = 26.0 



K^. = ( 1.00 - 0.60) g^gg ^ ^^a ^ g^gg ^ 

VO.86/' \v 1.32 J 



la = 0.333 + log/^iill^ = 0.761 

fron figure 2 . 
* 29.4 

and figure 5 gives 

\m = 2,42 X 10* 

Tor the lower surface; 

I 

, 5.5 (0.161) -Q-e , Q^QgQ ^ 
W 6,35 (2.71) 

- 2.58 loge [l.l8 (0,223 XlO**') 1.12 (6.35) lO' 

K.K « .(Q'Q.^ - . 9'.^^) — ^ 6.35 XIO^ = 72.3 X- 10® 

^ f 1 ^ 1 ^ 

\1,12> 'yLlS/ ■ 
I, = 0.00 + log (~^^^ = 0.053. = 25.0 

Kgft = (1.00 - 0 60) 6.35 X 10« = 8.06 X 10^ 

VO.86/ \1.18/ 
I- = 0.286 + log ^i—l^ = 0.602, Si » 28.8 
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and fron figure 5 * 
Xl » 1,92 X 10* 

The drag coefficient is therefore 

^ 3 (8.48 . 1.92) XIO* ^0.86)"* = 0.0098 
6,35 X 10® 

The profile drag of moat airfoils can be calculated 
ty the use of the operations involved in these two 
exeimples, 

DISCUSS lOR Airs CONCXUSIONS 



The graphical method for analyzing the turbulent 
region, presented and illustrated in the previous sec- 
tions, uses a rough approximation to the velocity dis- 
tribution. It is possible to use such an. appro ximat ion 
because the growth of the turbulent boundary layer seems 
rather insensitive to the actual shape of the velocity 
distribution in a region over which the velocity gradi- 
ent does not change markedly. This point can be 
illustrated by cons ider in^ the upper ^surface velocity 
distribution of the *IACA 35-315 airfoil in figure 6, It 
is seen that the tnib\ilent region has* three markedly 
different regions. 'Vh«n the cax'^n la^/ 1 :n is based on 
these natural divisions, the vaT..''o ^^ynd for the profile~ 
drag coefficient xs 0 0140 fexa-fi-Lt- l). Purther division 
of the velocity <2 istr""' bat ioi fa:j.t» to change this value. 
On the other hand, if the 'aacaral divisions are not 
taken into account, the cali-ulatad drag coefficient is 
considerably in error. For exajnple, if the above 
turbulent region is treated as a single segment, the 
drag coefficient is 0.0186, and when treated as two Beg- 
ments tho drag coefficient is 0.0159. 

The theoretical velocitjr distributions obtained by 
the method of reference 2 have the shortcoming of having 
a stagnation point at the trailing edge. In order to 
compensate for this incorrect trailing-edge velocity, it 
is necessary to. extrapolate the velocity gradient from 
ahead of about x/c " 0.'9, The e-rror involved in such 
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an approximation is rather small "because changes in the 
skin friction, due to trailing— edge velocity variations, 
are conpensated for Toy the corresponding changes in the 
pressure drag. This can he illustrated by assuming for 
the trailmg-edge Telocity in figure 6 a value of V/Vq 
= 0.90. The calculated profile— drag coefficient for this 
new trailing-edge velocity is 0.0136, which is in satis- 
factory agreeineiit with the previously calculated, value of 
0.0140. 

Che experimentally measured velocity distrihutions 
and profile-drag coefficients for a wide range of lift 
coefficients of the ffiOA, 35—315 airfoil are presented 
in reference 6. Theodorsen's method, reference 2, also 
was used to calculate the theoretical velocity distri— 
hut ions for the UAOA 35-215 airfoil at a numher of lift 
coefficients. The calculations were made at- specific 
lift coefficients rather than specific angles of attack 
because the distributions thus obtained more nearly aj)— 
proxinote experimental measurements. The profile— drag 
coefficients at these lifts were calculated from both 
the theoretical and the ea^jer imental velocity distribu- 
tions. She values thus obtained are compared with the 
experimentally measured polar in figure 8* The position 
of transition for the loiv— drag range of lift coefficients 
was found experimentally to lie from 5 to 10 percent of 
the chord length behind the position of maximum velocity. 
The calculated low^drag coefficients assume transition 
to occur at the position of maximum velocity and conse— 
q_uently these computed values are higher than the ex— 
perinental values. 

The -)rof lie— drag coefficients and experimental pres- 
sure distributions over a wide range of Kach numbers, 
and corresponding Seynolds numbers, are given m refer- 
ence 5 for an HAGA 66,2-420 airfoil. The procedure of 
the present report has been used to calculate the drag 
at varioxis I.'ach numbers for the smooth airfoil with tran- 
sition *ixed at x/c = 0.10, both at 0° angle of at- 
tack. These calculated values are compared with the 
experiment ally measured profile drags m figure 9. The 
values for the smooth' airfoil are in good agreement with 
experimental measurements up to about II = 0,53, which 
corresponds to a Reynolds number gf 16 X 10^. Above 
this Ilach number the measured drag increases rapidly be- 
cause transition moves ahead of the posit:i,on of maximum 
velocity* Since the turbulence level of the tunnel is 
higher than that of free air, the boundary— lay er Seynolds 
number at transition is well below the value of 8000 which 
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was mentioned as a criterion for transition on a smooth 
airfoil m flight, Moreover, as is shown m reference 7, 
compressiTjility effects produce a profound change in 
■boundrrjr-layer Eeynolda number at trans"ition. So, lack- 
ing Icnoirlodge of the proper transition criterion to use 
for the \7ind tunnel in which these tests were m^d-e, 
t''ansit:o.i on the smooth, airfoil was assumed to occur at 
the na::xnuzi velocity position. The calculated values 
for the case of transition fixed at x/c = 0,10 are 
unifornly lower than the experimental points "because the 
theory does not account for the drag of the cartorundum 
used to fix transition. When transition was fixed at 
x/c « 0.60, the experimentally measured drag coefficient 
was aljovit 0.0009 abov^ the measured values for the 
smooth airfoil. This drag increase is due chiefly to 
bhe crrhorundum strip. Since transition at x/c » 0.10 
was fixed Tjy a similar carhorundum strip, it is to "be 
expected that the drag of the carborundum in this case 
also uill be about 0,0009. Vhen this amount is added to 
the theoretically calculated drag, the calculated values 
are m satisfactory agreement vrith the experimental 
measuroaents. 

As a further comparison of the results predicted, 
by the .^ethods preseaiicsd here, with experimental measure- 
ments, ^Tolars for the iTACA 67,1—215 ai^-foil have been 
calculated for Reynolds numbers of 2,9 X 10^ and 
6.0 X 10^. The experimentally measured polars at these 
Reynolds numbers and a simple method for computing the 
ve?.ocity over low— drag airfoils at various lift coeffi— 
cionts ere presented m reference 8. In figure 10 the 
experinontal neasurcments are compared with the theoreti- 
cally calculated polars. The large discrepancies in the 
drag coefficients for high lift coefficients may be due 
to "ttitbulcnt separation" on the upper surface of the 
airfoil. Such a breakdown of the flow is believed to 
occur \rhon a thick turbulent boundary layer la subjected 
to a stoop ad'verse pressure gradient. To detormino the 
lijnits of applicability of the present method, polars 
were calculated for a number of widely different airfoils. 
Satisfactory agreement with experimental measurements was 
fovnd up to lift cooff ietents for which the drag coeffi— 
Gifent of the surface of highest drag on the airfoil was 
somewhac g^-eater than 0.01. With higher lift coefficients 
the ca?tcilated drag coefficients were progressively lower 
than the neasured values. It is also seen from figure 10 
that on incorrect variation of the lift coefficients with 
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Reynolds num'ber corresponding to the end points of the 
low— drag region is predicted ty the present method. 
Ihis shortcoming results from the lack of an adequate 
criterion for determining *he tranft it ion—point movement 
caused "by the velocity peak at the forward part of the 
airfoil, ^h is' peak ariees' from the superposition of 
small additional and haeic low-drag velocity distribu- 
tions. 

The problem of sudden shifts in transition point 
with snail changes in lift coefficient is not limited to 
low-drag airfoils. In figure 11 the calculated polar of 
the ISTACA 2S015 airfoil is compared with the experimental 
measurements given in reference 8. It is seen that the 
aTjrupt drag increase," m the wind tunnel at small negative 
lift coefficients occurs at a different lift coefficient 
from that predicted. Moreover, figure 11 indicates that the 
drag coefficients for conventional as well as low— drag 
airfoils can "be calculated for wide ranges of lift coef- 
ficients. 
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TABL3 I.- VALUES OP T ^ TOR OALCULATIITG t IN 

^0 KHy+S.lS 



A SSGION OP SECRMSIHO VBLOCIfY 





10« 


2x10® 


5X10* 


10^ 


2x10'' 


5x10'' 


"10"^" 


18.5 


0.0707 


0.0417 


0.0177 


0.0088 


0,0044^ 


0.0019 


0.0009 


19.5 


,1220 


.0666 


.0287 


.0145 


.0073 


.0030 


.0015 


80.5 


.1827 


.1039 


. 0459 


.0837 


.0120 


.0049 


.0025 


31.5 


,2638 


.1571 


. 0723 


.0380 


.0195 


.008 0 


.0040 


22.5 


.3623 


.2300 


.1112 


.0601 


.0314 


.0130 


.0066 


23.5 


.4793 


.3210 


.1661 


.0951 


.0497 


.0209 


.0106 


24.5 


.6104 


.4305 


.2396 


.1403 


.0774 


.0332 


.0170 


25.5 


.7518 


.6555 


'.3321 


.2048 


.1176 


.0520 


.0270 


26.5 


.9006 


,6985 


.4420 


.2881 


.1738 


.0798 


.0425 


27.5 


1,0544 


.8383 


.5676 


.3901 


.2479 


.1208 


.0659 


28.5 


1.2116 


.9900 


.7050~^ 


.5086 


.3409' 


,1774 


.1001 


29.5 


1,3708 


1.1458 


.8508 


.6406 


.'.520 


.2520 


.1485 


30.5 


1.5313 


1.3041 


1.0023 


.7824 


, S:-V76 


.3454 


.8139 


31.5 


1.6930 


1.4640 


1.1581 


.9314 


.'^147 


.4560 


,2975 


32.5 


1.8551 


1.6251 


1.3166 


1.0852 


.d303 


.5813 


.3993 


33.5 


2.0175 


1.7867 


1.4766 


1 . 2421 


1.0118 


.7183 


.5173 



lABLB II.- TAI.7BS OP /; = PGR OAC.CUT4ATI3T& 1 

v/q ::Rr-6.i3 



IN A ESGIOir OP IirCfiSAS INO rBLOCII?! 
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Figure 2.- Graph for evaluating turbulent 
^-factor in regions of decreas- 
ing velocity. 
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Figure 3.- Graph for evaluating turbulent ^-factor ^ 
m regions of constant velocity. 'j' 




Figure 4.- Graph for evaluating turbulent ^■ 
factor in regions of increasing 

velocxty. 
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Flgore 7.- Approximation to expeTinental velocity distribution for MJCk 66,2-420 airfoil at 
transition fixed at x/c s 0.10. 
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Pxgore 8.- Comparison cf calculated and experimentally measured polars for KaCA 35-215 airfoil at 
H = 5.75 X 10''. 
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ffigure 9,- CoiBpanson of calculated profile dra^ coefficients Tiith 

measured values at varioxis Mach numbers for HkCk 66,2-420, 

ot - 0°, 
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Figure 10,- Comparisor. cf celciilatecL and experimentally measured polars 
for NfaCA 67,1-215 airfoil. 
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J'igare 11.- Coiupariscn of calculated and experimentally measured polars for JikCk 23015 airfoil at 
H = 6.9 X 106. 
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